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Abstract In this letter we use graphene bubbles to study the Raman spec- 
trum of graphene under biaxial (e.g. isotropic) strain. Our Gruneisen pa- 
rameters are in excellent agreement with the theoretical values. Discrepancy 
in the previously reported values is attributed to the interaction of graphene 
with the substrate. Bilayer balloons (intentionally pressurized membranes) 
have been used to avoid the effect of the substrate and to study the depen- 
dence of strain on the inter-layer interactions. 
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Manuscript text Graphene is the latest carbon allotrope to be isolated^, 
and it is now at the center of a significant experimental and theoretical re- 
search effortP®. In particular, near-ballistic transport at room temperature 
and high carrier mobilities'^^ make it a potential material for nanoelectron- 
ics. 

Strain can be used to tailor the electronic properties of graphend 14 * 15 l 
This allows one to make an all-graphene circuit where all the elements are 



made of graphene with different amounts and types of strain^. Furthermore, 
certain configurations of strain are equivalent to a magnetic field 1 ™, which 
can be very higtW Thus, it is essential to be able to probe strained graphene 
and to distinguish between different types of strains. 

Elastic and inelastic light scattering are powerful tools for investigating 
graphene^^. Raman spectroscopy allows monitoring of doping, defects, 
disorder, chemical modifications and edges^^HZB All sp 2 bonded carbons 
show common features in their Raman spectra, the so-called G and D peaks, 
around 1580 and 1360 cm -128 . The G peak corresponds to the E2 S phonon at 
the Brillouin zone (BZ) center (T point). The D peak is due to the breathing 
modes of six-atom rings and requires a defect for its activation^. It comes 
from TO phonons around the BZ K point^ and it is activated by an intra- 
valley scattering process. The 2D peak is the second order of the D peak. 
This is a single peak in monolayer graphene, whereas it splits into four bands 
in bilayer graphene, reflecting the evolution of the band structure^. The 
Raman spectrum of graphene also shows significantly less intensive defect- 
activated peaks such as the D' peak, which lies at ~ 1620cm -1 . This is 
activated by an intra-valley process i.e. connecting two points belonging to 
the same cone around K (or K')^. The second order of the D' peak is 
called 2D' peak. Since 2D and 2D' peaks originate from a Raman scattering 
process where momentum conservation is obtained by the participation of 
two phonons with opposite wavevector (q and -q), they do not require the 
presence of defects. Thus, they are always visible in the Raman spectrum. 

Strain can be very efficiently studied by Raman Spectroscopy since this 
modifies the crystal phonon frequency, depending on the anharmonicity of the 
interatomic potentials of the atoms. The rate of this change is given by the 
Gruneisen parameter (7). Several experimental and theoretical works stud- 
ied uniaxial strain in graphene 122 ^, showing that this leads to: i) softening of 
the modes for tensile strain; ii) splitting of the G peak for increasing strain; 
iii) broadening or splitting of the 2D peak, depending on the strain direction. 
However, there is significant discrepancy among the reported results. The 
Gruneisen parameters are difficult to study in uniaxial strain because they re- 
quire the knowledge of the Poisson ratio, which depends on the substrate and 
on the degree of adhesion between graphene and the substrate itself-. Fur- 
thermore, it is difficult to calculate the D and 2D Gruneisen parameters since 
uniaxial strain moves the relative position of the Dirac cones. This affects 
the 2D peak phonon because this mode is activated by an inter-valley Raman 
scattering process^. Finally, the G and 2D Gruneisen parameters strongly 
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depend on the direction of the uniaxial strain and on its strength^®. 

Biaxial strain is more suitable for the calculation of the Gruneisen pa- 
rameter because it does not depend on the Poisson ratio and effects due to 
the relative anisotropic shifts of the Dirac cones are absent^. The Raman 
spectrum of graphene under biaxial strain has been experimentally studied 
only in two works: Ding et al.^ used a graphene flake deposited on a piezo- 
electric substrate. However, here the strain was not directly measured, but 
derived from the G peak position, using the calculated Gruneisen param- 
eter in Ref.^. The authors also found different Gruneisen parameters for 
the D and 2D peak, in contrast to the theoretical predictions^. Metzger et 
aP^ used a suspended graphene over a shallow depression. The strain was 
directly measured by Atomic Force Microscopy (AFM). In both works the 
Raman spectrum of graphene does not show any splitting of the G peak and 
the modes soften for increasing strain. However, very different Gruneisen 
parameters have been reported, Table 1. We will show that this discrepancy 
could be due to the strain or doping initially imposed on graphene by the 
substrate. 

In this letter we analyze the Raman spectrum of graphene under biaxial 
strain by studying graphene bubbles formed during the deposition of large 
graphene flakes on a oxidized silicon substrate (Si/SiO x ). Single layer and bi- 
layer bubbles have been measured and mapped by Raman spectroscopy and 
AFM. From the shift of the Raman peaks, we derived the Gruneisen parame- 
ters for G, D, 2D and 2D' peaks. Furthermore, we also studied the evolution 
of Raman spectra of bilayer graphene with biaxial strain by pressurizing a 
graphene membrane with nitrogen gas. Here bilayer is specially chosen to un- 
derstand the effect of biaxial strain on inter-layer interaction of AB- stacked 
bilayer graphene. Moreover this technique allowed us to eliminate any effect 
from strain or doping imposed on graphene by the substrate. 

Large single layer graphene flakes (lateral size above 0.1mm) were pro- 
duced by micro- mechanical cleavage of bulk graphite and deposited on Si/SiO a 
previously cleaned by oxygen plasma^. The bubbles can be easily identified 
by optical microscope because they are characterized by Newton rings, i.e. 
fringes of different colors, produced by the interference of the reflected and 
transmitted light waves between graphene and the silicon^. The bubbles 
have different shape and size^ we selected only large and spherical bubbles 
(average diameter of 5-10 /im). Graphene and bilayers bubbles have been 
identified by optical contrast^^ and by the shape of the 2D peak in the 
Raman spectrum^. 
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Figure 1: (color on-line) AFM image of a single-layer graphene bubble. The 
lateral scale is 4 fim. 

Atomic Force Microscopy has been used to measure the strain (e). Figure 
1 shows an AFM image of a bubble. The strain has been calculated in the 
following way: from the AFM we know the height (h) and width (w) of the 
bubble. Assuming that the strain is zero when there is no bubble (h=0), 
the strain will be: e = {L — Lq)/Lq, where L is the length of the arc after 
deformation and Lq is the initial length of the arc, i.e. Lq = w. We found that 
for most of the graphene bubbles e ~ 1%. From the shape of the bubble, it is 
possible to calculate the differential pressure by using the following relation: 
Ap = f(i>)Et[h 3 /(w/2) 4 ], where f(u) is a function of the Poisson ratio v, E 
is the Young's modulus and t is the thickness of the bubble. We used Et=347 
N/m and i/=0.16, thus f(u) = 3.09 40 . We estimated an average pressure of 
about 15 bar for most of the bubbles. 

A Witec Raman spectrometer, equipped with 488nm and 514. 5nm laser 
lines and a objective with numerical aperture of 0.95, has been used. This 
gives a laser spot size of about 400nm. The laser power is kept well below 
lmW. Careful control of the laser power is necessary: the bubble can move 
or disappear if the laser power is too high because mass transport under 
the bubble is activated by laser heating^. The Raman spectrometer is also 
equipped with a piezoelectric stage, which we used to Raman map the bubble 
and to take measurements exactly at the center of the bubble. The peaks have 
been fitted with a single Lorenzian lineshape and we analyzed the position 
and Full Width at Half Maximum (FWHM) of the peaks. 

The Gruneisen parameter has been calculated as^ 7 = [u — oo }/[2eu! }, 
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Figure 2: (color on-line) Schematic showing the experimental setup used to 
make the balloons; the inset shows an optical micrograph of a large bilayer 
covering a 50/im aperture in a copper scaffold. 

where u and uq are the Raman frequencies at finite strain and zero strain, 
respectively. 

Bilayer graphene balloons were produced by pressurizing a specially made 
metallic container covered with large bilayer graphene membranes with ni- 
trogen gas. Large bilayer graphene membranes were prepared by the method 
previously reported^. In brief, graphene crystals were exfoliated on top 
of PMMA coated (90 nm) silicon wafer by mechanical exfoliation. A thick 
copper film was deposited on top of the selected flakes by using a series 
of photolithography and electrodeposition steps. The deposited copper film 
contained an open aperture of 50/im in diameter, which was aligned with the 
chosen graphene crystal so that graphene fully covered the aperture. Cop- 
per scaffold with graphene layers were released from PMMA by dissolving in 
acetone. The samples were finally dried in a critical point dryer to prevent 
membrane from rupturing because of the surface tension. These graphene 
membranes were used to seal a small opening in a specially made metallic 
container with incorporated gas inlet and pressure gauge, Figure 2. Silver 
epoxy was used to attach the scaffold with graphene to the metallic container. 
Pressure inside the container was changed controllably by using nitrogen gas. 
We used a Renishaw Raman spectrometer (514. 5nm) to monitor the changes 
in the Raman spectrum of the balloon under various differential pressures 
(dp = Pinside — Pout)- The changes in Raman spectra of graphene balloons 
with different pressures were spatially uniform all over the membrane and 
were completely reversible (see Supporting Information). We were able to 
apply a maximum differential pressure of 2 bar. At higher pressures the 
membranes have a tendency to crack or burst. We speculate that this upper 
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limit of pressure was not limited by intrinsic property of graphene rather due 
to the weak adhesion of graphene to the metal scaffold or due to the extra 
strain at the edge of the aperture. 

Figure 3(a) shows the Raman spectrum of graphene measured on two 
adjacent regions: on the substrate, far from the bubble, and at the center of 
the bubble. A very small D peak is visible in both spectra, as shown in the 
inset. We can observe that the Raman spectrum taken at the center of the 
bubble is strongly blue-shifted compared to the Raman spectrum measured 
on the substrate. No splitting of the G peak and strong broadening/splitting 
of the 2D peak are observed on the bubble. In particular, the G peak position 
decreases from 1598 to 1525 cm -1 , while the 2D position moves from 2695 
cm -1 to 2552 cm -1 when moving from the substrate to the center of the 
bubble, Figure 3(b,c). The 2D FWHM stays constant at about 30 cm -1 . 
These changes in the Raman spectrum are in agreement with biaxial strain^. 
Note that the shifts in the peak position cannot be explained by doping, as 
the G peak position for undoped graphene is about 1582 cm -1 and it increases 
for increasing doping leveP^. 

We can now calculate the Gruneisen parameters by using the initial 
peak frequency measured for unstrained and undoped graphene. We used 
the following parameters (taken at 488nm): ojq(G)= 1582cm -1 ; Uo(2D) = 
2692cm" 1 ; u (2D') = 3245cm" 1 and u Q (D) = 1 349cm" 1 >™ We then find the 
following peak shifts: -57cm~ 1 /% for the G peak, -140 cm _1 /% for the 2D 
peak; -108 cm _1 /% for the 2D' peak and -68 cm _1 /% for the D peak. They 
give the following Gruneisen parameters: j(G)— 1.8; 7(2/})= 2.6; 7(2/}') = 
1.66 and ^y(D)= 2.52, in excellent agreement with previous calculations^, 
Table 1. In contrast to Ref.^ we do not observe any strong difference be- 
tween 7(2.D) and 7(-D) and no splitting of the D peak has been observed, 
inset Fig. 3(a). On the other side, our r y{2D') is in good agreement with the 
results reported in Ref.^. 

In the calculation of the Gruneisen parameters we did not use as initial 
peak frequencies the values measured on the substrate because it is well know 
that graphene on Si/SiO^ can be strongly doped^R Recent experiments have 
suggested the presence of a water layer between graphene and the substrate, 
which could cause strong p-doping^Sl. This layer should be absent under the 
bubble^. Indeed, the Raman spectrum measured at the center of the bubble 
shows a very small G peak intensity, compared to the 2D intensity, typical 
of undoped graphen d 18 * 44 l In contrast, the Raman spectrum of graphene 
on the substrate shows: a very high G peak position (1598 cm -1 ) and an 
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Figure 3: (color on-line) (a) Raman spectra measured at the center of a 
single-layer bubble and on the flat area on the substrate at 488nm;(b) Raman 
map of the G peak position; (c) Raman map of the 2D peak position. The 
2D and G peak positions strongly decrease when moving from the substrate 
to the center of the bubble. 
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Figure 4: (color on-line) a) Raman spectra measured at the center of the 
bilayer bubble and on the flat area on the substrate at 488nm; (b) Raman 
spectra of a bilayer graphene balloon measured for different differential pres- 
sure at 514. 5nm; (c) G and 2D peak position of a bilayer graphene balloon 
under different differential pressure. 8 



intensity ratio between G and 2D peak of about 1. This further confirms that 
graphene is strongly dopecP^ 44 1 45 1 . By using 1598 cm -1 as G peak position 
for unstrained graphene, we would get a Gruneisen parameter of 2.3, giving 
an error of 30%. 

Table 1 includes the most recent results on biaxial strain in graphene. We 
also report the ratio between the phonon peaks shift (A). This parameter 
does not require knowledge of the strain, so it can be used to compare the 
results of different works. 

We now focus on bilayer graphene bubbles. Figure 4(a) shows the Ra- 
man spectrum of such an object, measured on two adjacent regions: on the 
substrate, far from the bubble, and at the center of the bubble. No D peak 
is visible in both spectra. As seen for the graphene bubble, the Raman 
spectrum taken at the center of the bilayer graphene bubble is strongly blue- 
shifted compared to the Raman spectrum measured on the substrate. The 
G peak moves from 1590 to 1560 cm -1 . Note the shape of the 2D peak: the 
Raman spectrum measured on the bubble shows the typical features of an 
AB-stacked bilayer graphene^, in contrast to the Raman spectrum measured 
on the substrate. Here the 2D peak is almost a broad single-peak and the G 
peak intensity is very high, compared to the 2D intensity. This clearly shows 
that the bilayer graphene on the substrate is doped®'. Effect of strain im- 
posed by the substrate is also possible. Thus, in order to calculate the peak 
shifts, we used the following values: 1580 cm -1 for the G peak and 2664.5, 
2696, 2715 and 2728 cm -1 for the four 2D components^. We then found a 
shift of - 20cm _1 for the G peak and -50cm _1 for the 2D peak. This gives a 
shift ratio between 1.9 and 2.5, depending on the 2D phonon component. 

In order to avoid any effect on the Raman spectrum from doping or 
strain imposed by the coupling with the substrate, we investigated the Raman 
spectrum of bilayer graphene balloons. Figure 4(b) shows the first and second 
order Raman spectra measured at different pressure. Figure 4(b,c) shows that 
the Raman spectrum rigidly blue-shifts for increasing pressure; furthermore, 
the 2D peak shape does not change with the pressure. We found that the 
shift in the G peak and 2D peak is 34 cm -1 /bar and 73-79 cm -1 /bar for 
the 2D components, Figure 4(c). This gives a peak shift ratio between 2.15 
and 2.32, in good agreement with the results of the bilayer graphene bubble. 
Assuming a Gruneisen parameter of 1.8 for the G peak, we find that the 
maximum strain reached in the balloon is about 1.2%, which is too small 
to break the balloon or to strongly alter the electronic structure of bilayer 
graphene. This value is also in good agreement from the strain calculated 
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Table 1: Table reporting the main results obtained for biaxial strain in 
graphene. Note that Ref.^do not measure directly the Gruneisen parameter 
of the G peak, but use the value reported in Ref.^R The error bars take into 
account the uncertainty in the peak positions. A(2D/G) is the ratio between 
the shift of the 2D and G peaks. 
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by using w/2 = 25/zm and h=3.3 /im, where h is obtained from the formula 
linking the differential pressure with the bubble shape (here Et= 694 N/m 
being a bilayer 251 ). 

In conclusion, we have used graphene bubbles to study the Raman spec- 
trum of graphene under biaxial strain. We have calculated the Gruneisen 
parameters of the G, 2D, D and 2D' peaks: they are in excellent agreement 
with previously calculated parameters. We have also investigated the Ra- 
man spectrum of bilayer graphene bubbles and balloon: we found that in 
both cases the strain does not alter the AB-stacking configuration, probably 
because the strain is too small (~ 1.2%). We have shown that initial dop- 
ing or strain caused by the interaction of graphene with the substrate can 
strongly affect the measured Gruneisen parameters. This could explain the 
discrepancy in the previously reported Gruneisen parameters. 
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